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ABSTRACT 

We propose to measure the weak cosmic shear using the spatial derivatives of 
the galaxy surface brightness field. The measurement should be carried out in Fourier 
space, in which the point spread function (PSF) can be transformed to a desired form 
with multiplications, and the spatial derivatives can be easily measured. This method 
is mathematically well defined regardless of the galaxy morphology and the form of 
the PSF, and involves simple procedures of image processing. Furthermore, with high 
resolution galaxy images, this approach allows one to probe the shape distortions of 
galaxy substructures, which can potentially provide much more independent shear 
measurements than the ellipticities of the whole galaxy. We demonstrate the efficiency 
of this method using computer-generated mock galaxy images. 
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1 INTRODUCTION 

The coherent distortions of background galaxy im- 
ages by the intervening metric perturbations provide 
us a direct probe of t he large scale mass d is tribu- 
tion (see revie w s by Bartelmann fc Schneider! l200ll : 
IWittmanl 120021 : iRefregierl l2003h . Recently, several 
groups have claimed positive detections of the weak 
lensing effect and ob tained usefu l const ra ints on the 
cosm ol ogical model ([Bacon et al 120001: IKaiser et all 



200d: Ivan Waerbeke et all |2000|: 



Maoli et al 



2001 



Bacon et al 
— TT — 



2003; 



Wittman et al 



20001: 



Rhodes et all 2001 : van Waerbckeet all 



200 ll; iHoekstra et all 120021: iRefreeier et all 1200a 



Brown et al 

TT — 



200; 



2003 1; I Jarvis et alj|2003l; iRhodes et al.ll20o4 [Hevmans et al l 
2003 : IMassev et all 120051: lyan Waerbeke et all 120051 : iDahlei 
20061; iHetterscheidt et all l200fj: [Hoekstra et alj l200rj: 
Jarvis et al. 2006: ISchrabback et al.l l2006: Semboloni et all 
20061 ). In future weak lensing observations (e.g. VST- 
KIDS, DES^, VISTA darkCAM, Pan-STARRS3, LSST0, 
DUNE, SNAlfl, JDEJV0), if the photometric redshift can 
be well calibrated, we will be able to study the dark 
energy properties (its abundance and equation of state) 



Hamana ct al. 



E-mail:jzhang@astro. bcrkeley.edu 
see www.darkenergysurvey.org 
see pan-starrs.ifa.hawaii.edu 
see www.lsst.org 
see snap.lbl.gov 
see destiny.asu.edu 



using the redshift dependence of t he shear fields dHi 
12003 : lAbazaiian fc Dodelsonl 120031: Ijain fc Tavlori 12003' 



Bernstein fc Jai n 2004; H u fc Jainll20ol:ISong fc Knoxl2004 
Takada fc Jainl 12004 iTakada fc Whitel |2004|; llshakl booll 
Simpson & Bridle 2005; Zhang ct al. 2005; Hamicstad ct al. 
20061 : ISchimd et al.l l2006: Taylor ct al. 2006; Zhan 200$ . By 



constraining the growth factor of the mass perturbation and 
the geometrical distance as functions of redshift separately, 
weak lensin g provides a consistency check of the cosmologi- 
cal model teratochvil et all 12004 ISimpson fc Bridle! 120051 : 
IZhang et al.ll2005l : iKnox et all 12009 1 and opens a window 
for testing alternative gravity th eories l|Acquaviva et all 
120041 : ISoneil2005l : llshak et allboOfj ). 



An important and challenging job in weak lensing 
is to measure the weak cosmic shear (of order a few 
percent) from the shapes (or ellipticities) of the background 
galaxy images, which have large intrinsic variations. The 
existing methods are all based on convoluting the galaxy 
images with some weighting functions, an d are called 
the IN TEGRAL method s here aft er (see iTvson et all 
' l99d: iBonnet fc Mellieil 1 19951: IKaiser et all Il995l 



Luppino fc Kaiser 1997; Hoekstra et alJ 1998; Rhodes et all 



20001 : [Kaiserl l2000l: iBridle et all l200ll: IrBernstein fc Jarvisl 



20021: iRefregier fc Baconl 120031: IMassev fc Refregierl booa 
Kuiikenl 120061 : iNakaiima fc Bernstein! I2OO6I) . The INTE- 
GRAL methods typically have disadvantages in three 
aspects: 1. since the galaxy images are smeared by the PSF 
(either instrumental or environmental), the INTEGRAL 
methods involve at least two folds of convolutions, the math 
of which is complicated; 2. the details of the methods are 
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often sensitive to the galaxy morphology and the form of the 
PSF; 3. the shear information from the shape distortions of 
galaxy substructures is not considered. Strictly speaking, 
the shapelets method (see, e.g. ,Refregier 2003) may not 
be called an INTEGRAL method, because the galaxy 
weighting functions form a complete set of orthonormal 
shapelets which have very convenient mathematical prop- 
erties. It also has the potential of measuring the cosmic 
shears on galaxy substructures. However, since this method 
requires calibrations of the intrinsic distributions of the 
shapelet coefficients, it has strong dependence on the galaxy 
morphology. 

In this paper, we propose to use the spatial derivatives 
of the galaxy surface brightness field to measure the cosmic 
shear. This method was first used by Seljak and Zaldar- 
riaga (1999) on CMB lensing. We generalize their analysis 
by including the PSF and carrying out the measurement 
in Fourier space. This approach is well defined regardless 
of the galaxy morphology and the form of the PSF, and 
involves simple image processing procedures. Given a high 
image resolution, the method can potentially probe the cos- 
mic shear from galaxy substructures, greatly suppressing the 
shape noise. 

We begin by introducing the method in fJ21 In SJ31 
this approach is shown to work well on different types of 
computer-generated mock galaxy images with general forms 
of PSF. A brief summary is given in 311 



2 THE METHOD 

We derive the relation between the cosmic shear and the spa- 
tial derivatives of the galaxy surface brightness field without 
a PSF in £|2 . 1 1 In the presence of an isotropic Gaussian PSF, 
the relation is modified and shown in i]2.2l In £|2.3I Fourier 
transformation is introduced not only to simplify the mea- 
surement of the spatial derivatives, but also to deal with 
general forms of PSF. 

2.1 Without the PSF 

The surface brightness on the image plane fi{0 ! ) and on the 
source plane fs(9" S ) [9 1 and 9 s are the position angles on 
the image and source plane respectively) are related through 
a simple relation: 

fi(o I ) = fs(e s ) 



(i) 



A0- 



where Ay = 5y + <&y, and $y = dS8( /dOf are the spatial 
derivatives of the lensing deflection angle, which can be ex- 
pressed in terms of the convergence k — ($ mx + & yv )/2 and 
the two shear components 71 = (<& xx — <& yy )/2 and 72 = 
Using eq.[T], we get: 



oej de' aef 



(2) 



dfs 



Of 



where we have implicitly assumed that <E>y is small, which 
is true for weak lensing. Assuming the original surface 
brightness field fs is isotropic on the source plane, the 



quadratic combinations of the derivatives of the lensed 
image provide a direct me asure of the shear components 
(|Seliak fc Zaldarriagalll999l ): 



1 {{d x fxf - (dyfjf) 

2 <(aji) 2 + {dyfry) 

{dxfidyfi) 



= -71 



(3) 



where the averages are taken over the whole galaxy. 

2.2 With an Isotropic Gaussian PSF 

The presence of PSF brings both advantages and disadvan- 
tages. On the positive side, the PSF smooths out the galaxy 
surface brightness field, which is originally not differentiable 
due to structures on arbitrarily small scales. On the other 
hand, the convolution of the galaxy image with the PSF 
leads to a nontrivial modification to eq. [3J, the form of which 
is calculated in this section. For simplicity, we assume the 
PSF is isotropic and Gaussian. General forms of PSF will 
be discussed in N2.3I 

The observed galaxy surface brightness distribution fo 
is related to // via: 



fo(9)= / d 2 9 1 w (9~9 I )f I (9 I ) 



where Wp is the Gaussian PSF with scale length j3: 



Wp{0) = 



2tt/3 : 



■ exp 



2/3 2 



(4) 



(5) 



Using eq.p] to replace // with fs and 6 1 with 8 s in eq. [4], 
we get: 



fo{e) = |det(A)| / d 2 e s W P {6- A6 s )f s (e s ) 



or equivalently: 



(6) 



fo{A0) = |det(A)| / d 2 9 s W f 3[A(8-8 s )]f s (9 s ) (7) 



= |det(A)| J d 2 ^ f s (9 s )W (8-9 s ) 

x [l-(9-9 s )-(A-I)-{9-9 s )//3 2 ] 

where I is the 2x2 unitary matrix. Note that the second 
part of eq.[7] is a result of Taylor expansion of the term 
W/3[A(9 ~ 9 s )] due to the small amplitudes of the lensing 
components <J?y. For convenience, let us define: 



Fs(9)= I d 2 9 s f s {9 s )Wf30 - 9 s ) 



(8) 



which is the surface brightness field we would observe in 
absence of lensing. Eq.[7] can then be re- written as: 



|det(A)| 



= (I - $xx - $yy)Fs(9) 



3 2 / „ d 2 F s , nA d 2 F s 



(9) 



89 2 



+ 2$ xy 



09 X d9y 



38% 



Let 8° = A^, then: 
dfo(9°) _ (A -i ) . dfo(A9) 



08? 



(10) 
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Using eq.[9] and eq.JTO], it is not hard to express the deriva- 
tives of fo in terms of the derivatives of Fs: 

dxfo 



|det(A) 

dyfo 
|det(A) 

where 
F ■ ■ 



= (1 - 2$ X X - <byy)F x - $ X yF y 

/3 (J&xxFxxx H~ 2-*&xy Fxxy H~ ^yyFx 



(ll) 



= (1 - 



$ i.^yy -^yyy ^^xy-^xyy ~\~ ^xi-fxxy) 

d n F s 



(12) 



Note that we have implicitly assumed that the spatial fluc- 
tuation of the cosmic shear is negligible on galactic scales. 
Assuming the distribution of Fs is isotropic, we obtain the 
following relation between the shear components and the 
spatial derivatives of the surface brightness field: 



1 ((d*fo) 2 - (d y fof) 

2 ((d x fo) 2 + (dyfo) 2 ) + A 

(d x fodyfo) 

((d x f y + (dyfo) 2 ) + A 

where 

A= y<V/o-V(V 2 /o)> 



= -7i 



= -72 



(13) 



(14) 



The derivation of eq. [13] and eq.[l4] is shown in the 
appendix. Note that in the limit when the galaxy image 
is very smooth over the scale length /?, the correction A 
approaches zero, eq.[T3] then reduces to eq. [3]. 



2.3 Fourier Transform and General PSF 

For the method to become useful, there are at least two 
remaining issues to be addressed: 1. how to measure the 
spatial derivatives of the surface brightness field; 2. how to 
deal with other forms of the PSF. It turns out that Fourier 
transformation provides a solution to both problems. 

Since convolutions in real space correspond to multipli- 
cations in Fourier space, one can easily transform the PSF 
to a desired form (an isotropic Gaussian form in our case) 
by multiplying the Fourier modes of the observed image 
with the ratios between the Fourier modes of the desired 
PSF and those of the original PSF (known from calibra- 
tions with stars). This operation is usually well defined if 
the scale length of the desired PSF is larger than that of the 
original PSF. Moreover, it turns out that for the purpose of 
measuring the cosmic shear, one does not need to transform 
the new image back to real space, because the derivatives of 
the surface brightness field can be more easily measured in 
Fourier space. 

As an example, we show how to measure quantities such 
as (jV/| 2 ), where / is the surface brightness field of interest. 
First of all, the distribution / in real space should be sam- 
pled with an interval AO which is a few times less than the 
size of the PSF to avoid translating high frequency power 
into the frequency range determined by t he sampling resolu - 
tion through discrete Fourier transform (|Press et al.lll992T ). 
In other words, the galaxy image should be "oversampled" to 
avoid aliasing power from small scales in the discrete Fourier 



transform. For undersampled images, one can smooth the 
images with an additional large enough PSF, which can be 
treated as a part of the PSF from the instrumentation, and 
therefore does not affect our discussion below. Similarly, to 
avoid such aliasing power at low frequency, the box size for 
the Fourier transform should be a few times larger than the 
image size. Given this setup, the Fourier transform of the 
image is defined as: 



JV-l JV-l 



f(li,lj) = A6 2 ^2 ^2f(0-m,9n)exp[i(e m k+9nlj)] (15) 

m— n— 

where 



0m(n) = m(n) x AO, m(n) = 0, 1, N - 1 
l iU) = i(j) x Al, i(j) = -N/2,..,N/2 
Al = 2w/(NA0) 



(16) 



N is the box size, chosen to be a power of 2 for the Fast 
Fourier Transform. It is now straightforward to show that 
(|V/| 2 ) can be expressed as the sum over the Fourier modes 
weighted by the wave numbers: 



(17) 



1 



JV/2 JV/2 



jV 2 A6» 4 



\hh,h)\ 2 (i 2 +i. 



i=-N/2 j = -N/2 



Ed. [17] gives exactly the quantity (|V/| 2 ) multiplied by the 
number of bright pixels covered by the galaxy, because the 
dark pixels have no contribution^]. Similarly, we can calcu- 
late the other terms in eq. [13] in Fourier space. Note that 
for the purpose of obtaining 71 and 72, it is not necessary 
to calculate the number of bright pixels because it appears 
in both the nominator and the denominator in eq.|13|. 



3 THE TEST 

This section is organized as follows: in i]3.1l we test the 
method using mock regular galaxies smeared by different 
forms of PSF; in £13.21 using mock irregular galaxies gen- 
erated by 2-D random walks, we further demonstrate the 
usefulness of this approach on galaxies with a different mor- 
phology, and explore the possibility of suppressing the shape 
noise in the shear measurements by including the informa- 
tion from galaxy substructures. 



3.1 With Mock Regular Galaxies 

Each regular galaxy in our simulation contains a 
thin circular disk with an exponential profile and 
a co-axial de Vaucouleu rs-type spheroidal component 
|de Vaucouleurs et al.|[l99lh . When viewed face-on, the sur- 
face brightness distribution (before lensing and smearing by 
the PSF) of the galaxy can be parameterized as: 



6 In the presence of noise, extra procedures may be required to 
clean the galaxy map before the Fourier transform. We shall dis- 
cuss this in a future paper. 
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Figure 1. Two PSF models used in the simulations of 93,11 The 
contours mark 3.2%, 6.3%, 12.5%, 25%, and 50% of the peak 
intensity. The left and right one correspond to Wr and Wr 
(rotated by certain angles) in eq.|19| respectively. 



f(r) = exp(-r/r d ) + f a /dexp 



-(r/r s 



(18) 



where r is the distance to the galaxy center, r d (r s ) is the 
scale length of the disk (spheroid), and f s /d determines the 
relative importance of the spheroid. The overall luminosity 
of the galaxy is only important in the presence of noise, 
which will be discussed in a future paper. 

Our simulation box is 128x128. We choose Vd to be 1/32 
of the box size of the simulation, r s = r^/2, and f s / d — 1. 
Note that changing these particular numbers does not af- 
fect our main conclusions. Once the galaxy's face-on image 
is generated, it is projected onto the source plane with a 
random inclination angle along a random direction perpen- 
dicular to the line of 

The projected galaxy image is subsequently distorted 
by a constant cosmic shear and smeared by the PSF in real 
space. We consider two PSF models given by the following 
forms rotated by certain angles (shown in FigfT}: 

WP{x,y) « exp [-(\x -y\ + \x + y\f/(8r 2 )] (19) 
W^ 2) (x,y) ocexp [-{x 2 + 0.8y 2 )/(2r 2 )] 

where r is the scale length, which is equal to six times 
the grid size, comparable to the galaxy size. The shear 
components (71,72) are chosen to be (—0.012,0.035), 
(-0.032, -0.005), (0.01, 0.02) for W r (1) , and (0.015, -0.024), 
(0.05,0.01), (-0.04,-0.04) for W r (2) . 

To measure the cosmic shear, we follow the procedures 
described in H2.3I The desired PSF has an isotropic Gaus- 
sian form with a scale length about 4/3 times that of the 
original PSF. The results are plotted in Fig[2] The results 
are consistent within lcr error regardless of the form of the 
PSF. 



3.2 With Mock Irregular Galaxies 

Our irregular galaxies are generated using 2-D random 
walks. The random walk starts from the center of the simu- 
lation box for 20000 steps, each of which is equal to the grid 



7 The intrinsic flattening parameter q of the spheroid part is set 
to one for simplicity. 




-0.04-0.02 0.02 0.04 
7i 

Figure 2. The measured values of (71,72) (averaged over 10000 
mock regular galaxies) are plotted as red (from galaxies smoothed 
by W^ 1 ') and blue (from galaxies smoothed by W^) dots with 
(small) lcr error bars. The centers of the black crosses are the 
input values of (71 ,72). Note that to be general, and in 

the simulations have been rotated by certain angles with respect 
to their definitions in eg. 1191 . 



size of the simulation box (which is now 1024 x 1024). Once 
the distance from the center is more than 1/6 of the box 
size, the walk starts from the center again to finish the rest 
of the steps. The surface brightness of the galaxy is equal 
to the density of the trajectories. Note that these galaxies 
naturally have abundant substructures, which are useful not 
only for further testing the method, but also for illustrating 
how much lensing information may be contained in the sub- 
structures. We caution that our random-walk-type galaxies 
are not based on any physical models, therefore they do not 
necessarily mimic observed irregular galaxies. In a future pa- 
per, more realistic galaxy models will be adopted to study 
this topic. 

For the purpose of this section, we smooth the galaxies 
directly with the isotropic Gaussian PSF of different scale 
lengths, which correspond to different angular resolutions. 
The scale length (5 (defined in eq.[5]) is chosen to be 1/256, 
1/128, 1/64, and 1/32 of the box size (roughly correspond- 
ing to 1/85, 1/43, 1/21, and 1/10 of the galaxy size). Fig|3] 
shows typical images of our irregular galaxy under these 
four different angular resolutions. For convenience, we plot 
the minimum j3 as unity in the figures of this section. 

The shear component 72 is set to zero, and 71 is fixed 
at 0.03. After averaging over 10000 irregular galaxies, we 
find that the measured 71 is 0.0324 ± 0.0029 for /3 = 8, 
0.0301 ± 0.0022 for /3 = 4, 0.0291 ± 0.0015 for (3 = 2, and 
0.0293 ± 0.0010 for ,9=1. More interestingly, as shown in 
FigHl the statistical error bar is found to decrease signifi- 
cantly when the angular resolution is increased. This is fur- 
ther illustrated in Fig[5] which shows an approximate power- 
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Figure 3. Sample images showing our mock irregular galaxies 
generated using 2-D random walks under four different angular 
resolutions. 



law relation between the measured variance of 71 and /3, the 
exponent of which is close to one. Note that as the angular 
resolution increases, one gets additional information on the 
cosmic shears from the galaxy substructures. If we naively 
assume that each bright pixel on the galaxy map provides 
an independent measurement of the cosmic shear, we ex- 
pect the variance of the measured cosmic shear to scale as 
the inverse of the number of the bright pixels, or P d , where 
d is t he Hausdorff dimension of the galaxy image l|Hausdorfj 
Since the Hausdorff dimension of our random-walk- 
generated irregular galaxies is 2 l|Falconerlll986h , the vari- 
ance of 71 should scale as /3 2 , which is not too far from 
what we have observed in our numerical experiments. In re- 
ality, substructures generated by the 2-D random walks are 
correlated at some unknown level, therefore, the observed 
exponent indicated in FigO is less than the Hausdorff di- 
mension. 



4 SUMMARY 

We have presented a simple approach of measuring the weak 
cosmic shear using the spatial derivatives of the galaxy sur- 
face brightness field. The measurement should be carried out 
in Fourier space, in which it is easy to evaluate the spatial 
derivatives and to transform the PSF to a desired form. The 
accuracy of the method is demonstrated using computer- 
generated mock regular and irregular galaxies. We find no 
systematic errors on the measured shear components in the 
numerical experiments. 

Given high image resolutions, this new method may re- 
duce the shape noise in the shear measurement significantly, 
because it takes into account the shape information on the 
galaxy substructures. Using the mock irregular galaxies gen- 



Figure 4. The probability distribution functions (PDF) of the 
measured 71 for four different image resolutions. The plots are all 
generated using 10000 mock irregular galaxies. The peak of the 
PDF is fixed at 1 in each plot. 
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Figure 5. The variance of 71 for four different values of (3. 

erated by 2-D random walks, we have shown that the vari- 
ance of the measured shear is indeed suppressed by a large 
factor when the image resolution is increased. This exam- 
ple encourages us to test this method on real galaxies of 
a wide range of morphology classes in a future paper by 
joining the Sh ear TEsting programme (Hey mans et al . 2006; 
iMassev et al.l 120061 ) . the results of which may be useful for 
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optimizing the signal to noise ratio in shear measurements 
and planning future weak lensing survey. 
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APPENDIX - RELATING THE COSMIC 
SHEARS WITH THE SPATIAL DERIVATIVES 
OF THE SURFACE BRIGHTNESS FIELD 

From eq. [IT], we have: 

IdeW ~ i9yf0)2] m 

= ( 1 - 6k) (F 2 - F 2 ) - 2 7 i (F? + Fy ) 

- P 2 \2kUi +7i(A + Ti) + 72 (T 2 - A)] 



and 



d x fodyfo (21) 



I dot ( a; 

= 2(1 — Qk)F x F v — 2y 2 {F 2 + Fy) 
- P 2 [2 K n 2 + 7l (A + T 2 )+ 72 (A-T 1 )] 
where 

A. F x F X XX H~ F x F X yy ~f" Fy F X xy ~\~ Fy Fyyy (^^) 

A FyFxX'X ~\~ FyFxyy F X F X Xy F X Fyyy 

n 1 F X F X X X H~~ Fx F X yy Fy F X Xy Fy Fy y y 

-Tl-2 F X F XX y H - F X Fyyy -\~ FyFxyy "t - FyFxXX 

^1 F X F X XX 3F X F X yy 3FyF XX y ~\~ Fy Fyyy 

^2 FyFxXX ^FyFxyy ~\~ 3F X F XX y ~~ F X Fyyy 

Note that according to the definitions in eq.[22], A is a 
scalar, A is a pseudo scalar, III + iU.2 is a spin-2 field, and 
Ti + 2X2 is a spin-4 field. If the intrinsic surface brightness 
distribution is isotropic, the spatial averages of A, TTi , II2, 
Ti, and T2 must vanish. As a result of this, we have: 

\{{d*fof - (d v fof) (23) 
= -7i ((f! + F 2 } + ^ (A) 



(d x fodyfo) 

-., ( (F 2 + F 2 ) + t {K 



We have neglected the factor |det(A) | which is equal to unity 
to the 0th order. Using the fact that A = VFs ■ V(V 2 F S ), 
and Fs = fo to the 0th order, it is now straightforward to 
prove ea. [13]. 



© 2006 RAS, MNRAS 000, Qj47] 



